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Abstract Two aspects of quark matter at high density are addressed: one is color super- 
conductivity and the other is ferromagnetism. We are mainly concerned with 
the latter and its relation to color superconductivity, which we call color mag- 
netic superconductivity. The relation of ferromagnetism and chiral symmetry 
restoration is also discussed. 

1. Introduction 

Nowadays it is widely accepted that there should be realized various phases 
of QCD in temperature (T) - density (ps) plane. When we emphasize the 
low T and high pB region, the subjects are sometimes called physics of high- 
density QCD. The main purposes in this field should be to figure out the prop- 
erties of phase transitions and new phases, and to extract their symmetry break- 
ing pattern and low-energy excitation modes there on the basis of QCD. On the 
other hand, these studies have phenomenological implications on relativistic 
heavy-ion collisions and compact stars like neutron stars or quark stars. 

In this talk we'd like to address magnetic properties of quark matter at low 
temperature. We first discuss the ferromagnetic phase transition and then a pos- 
sibility of the coexistence of ferromagnetism (FM) and color superconductivity 
(CSC). We also present an idea about how FM is related to chiral symmetry. 
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CSC should be very popular and many people believe that it is robust due to 
the Cooper instability even for small attractive quark-quark interaction in color 
3 channel [1]. On the contrary, we are afraid that FM has not been so familiar 
yet. So, we'd like to begin with a brief introduction about our motivation for 
the study of FM. 

Phenomenologically the concept of magnetism should be directly related 
to the origin of strong magnetic field observed in compact stars [2]; e.g., it 
amounts to O(IO^^G) at the surface of radio pulsars. Recently a new class of 
pulsars called magnetars has been discovered with super strong magnetic field, 
Bg ~ 10^^~^^G, estimated from the P — P curve [3, 4]. First observations 
are indirect evidences for super strong magnetic field, but discoveries of some 
absorption lines stemming from the cyclotron frequency of protons have been 
currently reported [5]. 

The origin of the strong magnetic field has been a long standing problem 
since the first pulsar was discovered [2]. A naive working hypothesis is the 
conservation of the magnetic flux and its squeezing during the evolution from 
a main-sequence progenitor star to a compact star, Bg cx R^^ with R being 
the radius [6]. 

Table 1. Surface magnetic field and the radius of stars by the conservation of the magnetic 
flux. 





Bs[G] 


7? [cm] 


Sun (obs.) 


10^ 


IQio 


Neutron star 


10" 


10^ 


Magnetar 


IQis 


10^ 



The relation of the radius and the expected strength of the magnetic field is 
listed by the use of the hypothesis in Table. 1 . Then, it looks to work well for 

explaining the strength of the magnetic field observed for radio pulsars. How- 
ever, it does not work for magnetars; considering the Schwatzschild radius, 

Rsch = 2GM/c^ = 4 X 10^ [cm] » 10^ [cm], (1) 

for the canonical mass of M = 1.4M0, we are immediately led to a contradic- 
tion. 

Since there should be developed hadronic matter inside compact stars, it 
would be reasonable to consider a microscopic origin of such strong magnetic 
field: ferromagnetism or spin polarization is one of the candidates to explain 
it. To graspe a rough idea about how hadronic matter can give such a super 
strong magnetic field rather easily, it should be interesting to compare typical 
energy scales in some systems (see Table 2): the magnetic interaction energy 
is estimated as -Emag = Mi-^ with the magnetic moment, /Hj = ei/{2mi). Thus 
we can see E^p <^ -Emag for electrons, while E^yp > E'mag for nucleons 
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Table 2. Magnetic interaction energies -Emag for 10 G and the typical energy scales £typ in 
electron, nucleon and quarlc systems. 





electron 


proton 


quark 


m,[MeV] 


0.5 


10^ 


1- 100 


S^ag[MeV] 


5-6 


2.5 X 10"^ 


2.5 X 10"^ - 2.5 


Etyp 


KeV 


MeV 


MeV 



or quarks. This simple consideration may imply that strong interaction gives 
a feasible origin for the strong magnetic field. The possibihty of ferromag- 
netism in nuclear matter has been elaborately studied when the pulsars were 
observed, but negative results have been reported so far [7]. Here we consider 
its possibility in quark matter from a different point of view [8]. 

2. What is ferromagnetism in quark matter? 

Quark matter bears some resemblance to electron gas interacting with the 
Coulomb potential; the gluon exchange interaction in QCD is similar to the 
electromagnetic interaction in QED and color neutrality of quark matter corre- 
sponds to charge neutrality of electron gas under the background of positively 
charged ions. It was Bloch who first suggested a mechanism leading to fer- 
romagnetism of itinerant electrons [9, 10]. The mechanism is very simple but 
largely reflects the Fermion nature of electrons. Since there works no direct 
interaction between electrons as a whole, the Fock exchange interaction gives 
a leading contribution; it can be represented as 

l^Fock = -e -Tj r^, (2) 

|k — 

between two elections with momenta, k and q, and spin polarizations, C and 
C', where the vector ^ specifies the definite spin polarized state, e.g. C = 
(0, 0, ±1) for spin up and down state. Then it is immediately conceivable that 
a most attiactive channel is the parallel spin pair, whereas electrons with oppo- 
site polarizations gives null contribution. This is nothing but a consequence of 
the Pauli exclusion principle: electrons with the same spin polarization cannot 
closely approach to each other, which effectively avoid die Coulomb repulsion. 
On the other hand a polarized state should give a larger kinetic energy by rear- 
ranging the two Fermi spheres. Thus there is a trade-off between kinetic and 
interaction energies, which leads to a spontaneous spin polarization (SSP) or 
FM at some density ^ . One of the essential points we learned here is that we 
need no spin-dependent interaction at the original Lagrangian to see SSP. We 
can see a similar phenomenon in dealing with nuclear matter within the rela- 
tivistic mean-field theory, where the Fock interaction can be extracted by way 
of the Fierz transformation from the original Lagrangian [11]. 
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Then it might be natural to ask how about in QCD. We list here some fea- 
tures of QCD related to this subject. (1) the quark-gluon interaction in QCD is 
rather simple, compared with the nuclear force; it is a gauge interaction like in 
QED. (2) quark matter should be a color neutral system and only the exchange 
interaction is relevant like in the electron system. (3) there is an additional fla- 
vor degree of freedom in quark matter; gluon exchange never change flavor 
but it comes in through the generalized Pauli principle. (4) quarks should be 
treated relativistically, different from the electron system. 

The last feature requires a new definition and formulation of SSP or FM in 
relativistic systems since"spin" is no more a good quantum number in relativis- 
tic theories; spin couples with momentum and its direction changes during the 
motion. It is well known that the Pauli-Lubanski vector is the four vector 
to represent the spin degree of freedom in a covariant form, 



In the rest frame, 



W = -jwiV". (3) 



^. = 0. :^ = l,V7 = ^i:. (4) 



where S = ^ ^ ^ ^ in the usual basis. For any space-hke four vector a 
orthogonal to k, a^k^ = 0, we then have a property, 

W^-« = -^75#. (5) 
By taking a 4-pseudovector a'^ s.t. 

with the axial vector C„ we can see the operator 

P(a) = ^(l + 75{^) (7) 

is the projection operator for the definite spin-polarized states; actually a'^ is 
reduced to a three vector (0, C,) in the rest frame and we can allocate C, = 
(0, 0, ±1) to spin "up" and "down" states. Thus we can still use C, to specify 
the two intrinsic polarized states even in the general Lorentz frame. 

We briefly present a heuristic argument how quark matter becomes ferro- 
magnetic by the use of above definition [8]. The Fock exchange interaction, 
/kCqC' between two quarks is defined by 

X /ox 

/kCqC = TT-ET-A^kCqC- («) 
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A^kCqC is usual Lorentz invariant matrix element and can be written in the 
lowest order as 



where the spin dependent term renders 
1 



a - h = t; 



ml 



-(k-C)(q-0 + "^'C-C' 



+ {m(Efc + m)(C • q)(C' • q) + m{E^ + m)(C' • k)(C • k) 

+ (k-q)(C-k)(C'-q)}/(£;fe + m)(£:g + m)] . (10) 

It exhibits a complicated spin-dependent structure arising from the Dirac four 
spinor, while it is reduced to a simple form, 



2 .^ i + C- C 

9^ (k-q; 



in the non-relativistic Umit as in the electron system. Eq. (11) clearly shows 
why parallel spin pairs are favored, while we cannot see it clearly in the rel- 
ativistic expression (10). We have explicitly demonstrated that the ferromag- 
netic phase should be realized at relatively low density region [8]. 

Relativistic ferromagnetism 

If we understand FM or magnetic properties of quark matter more deeply, 
we must proceeds to a self-consistent approach, like Hartree-Fock theory, be- 
yond the previous perturbative argument. In ref. [11] we have described how 
the axial- vector mean field (AV) and the tensor one appear as a consequence 
of the Fierz transformation within the relativistic mean-field theory for nuclear 
matter, which is one of the nonperturbative frameworks in many-body theo- 
ries and corresponds to the Hatree-Fock approximation. We also demonstrated 
they are responsible to ferromagnetism of nuclear matter. An important point 
obtained there is the "condensation" of AV. ^ 

When we consider the non- vanishing AV in quark matter, 

V = -7573UA, Va//z, (12) 
we see an interaction between quarks and AV, 

Hint^CT-'^A = (T^UA, Ua>^, (13) 

in a similar form to the magnetic interaction in QED. Then the quark propaga- 
tor in AV renders 

G-y,\p)=i)-m-^ + -1^pA. (14) 
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The poles of Ga{p), detG^ (po=en)=0, give the single-particle energy spec- 
trum: 

en = ±e± (15) 

e± = ^p^ + \J\ + m^± 2^m?\]\ + (p • Ua)2, (16) 

where the subscript ± in the energy spectrum represents spin degrees of free- 
dom, and the dissolution of the degeneracy corresponds to the exchange split- 
ting of different "spin" states [10]. Actually it is reduced to a familiar form, 

2 

e-t = m + ^ ± Ua, in the non-relativistic limit, while 

e± = + {\pz\ ± Ua? (17) 

in the extremely relativistic limit, m ^ 0. Note that Ua only shifts the value 
of momentum in Eq. (17), so that it should be redundant in the massless case. 

There are two Fermi seas for a given quark number with different volumes 
due to the exchange splitting in the energy spectrum. The appearance of the 
rotation symmetry breaking term, oc p • Ua in the energy spectrum (16) im- 
plies deformation of the Fermi sea: so rotation symmetry is violated in the 
momentum space as well as the coordinate space, 0(3) 0(2). Accordingly 
the Fermi sea of majority quarks exhibits a prolate shape (F~), while that of 
minority quarks an oblate shape (F+) as seen Fig. 1 ^ . 






Figure 1. Modification of the Fermi sea as Ua is increased from left to right. The larger 
Fermi sea (F^) takes a prolate shape, while the smaller one (F^) an oblate shape for a given 
Ua- In the large Ua limit (completely polarized case), disappears as in the right panel. 



The mean spin polarization is then given by 



trbsTsOA^)] 



C (27r)4 

d^p Ua + P 
F+ (27r)3 e+ 



+ 



(27r)3 e_ 



(18) 



(19) 
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with (3 = \Jpj. + m^, from which we can immediately see the non-vanishing 
value of Ua gives rise to spin polarization. Since the spin polarization is not 
necessarily measurable quantity, we'd better to see another observable, the 
magnetization, which is defined as the magnetic moment per unit volume and 
the magnetic field directly couples with it. In QED the magnetic field couples 
with quarks by way of the term, fJ-qQcr^i^qF'^'^ , with the Dirac magnetic moment 
fiq = e/(2m), and we can easily see that the magnetization M is directed to 
the z direction; 



—I 

1 
2 



c (2^ 



tr[7o7573GA(p)] 



f (Pp m r 



(fp 



m 



F+ (27r)3 /? Jf- (27r)3 (3 



(20) 
(21) 



with P = ^Jp^^^rmr in the units of the Dirac magnetic moment. Note that 
the magnetization of each Fermi sea has now the opposite direction and it does 
not explicitly depend oxvUa, but the net magnetization arises by way of the 
exchange splitting of the Fermi sea. Thus we see the ground state holds ferro- 
magnetism in the presence of Ua- 



3. Color magnetic superconductivity 

If FM is realized in quark matter, it might be in the CSC phase. In this 
section we discuss a possibility of the coexistence of FM and CSC, which we 
call Color magnetic superconductivity [13]. 

In passing, it would be worth mentioning the corresponding situation in 
condensed matter physics. Magnetism and superconductivity (SC) have been 
two major concepts in condensed matter physics and their interplay has been 
repeatedly discussed [14]. Very recently some materials have been observed to 
exhibit the coexistence phase of FM and SC, which properties have not been 
fully understood yet; itinerant electrons are responsible to both phenomena 
in these materials and one of the important features is both phases cease at 
the same critical pressure [15]. In our case we shall see somewhat different 
features, but the similar aspects as well. 

We begin with an OGE-type action: 



-9 



D^u{x,y) 



,(22) 



where D'*'' denotes the gluon propagator. By way of the mean-field approxi- 
mation, we have 



Imf 



2j J2^\ 



(23) 
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in the Nambu-Gorkov formalism. The inverse quark Green function G^^{p) 
involves various self-energy (mean-field) terms, of which we only keep the 
color singlet particle-hole V{p) and color 3 particle-particle (A) mean-fields; 
the former is responsible to ferromagnetism, while the latter to superconduc- 
tivity, 



G-\p) 



- m + ^ + V{p) 7oA'l"(p)7o_ 
A(p) ;^ — m — ^ + V{p) 

Gnip) Guip) 

G2l{p) G22{p) J 

where 

^^{k) = CiP^{-k), V = CV^C-\ (25) 

Taking into account the lowest diagram, we can then write down the self- 
consistent equations for the mean-fields, V and A: 

-V{k) = {-igf j ^^{-iD^'^{k-p)}^,^{-^G^,{p)}^,^. (26) 



and 



= {-^gf / ^^{-^D^^[k-p)}^,^^{-^G2M}lu^■ 

(27) 



G2r 



Gi. 



Nc+l 
Nc 



Nc-1 



Nc 



Figure 2. Graphical interpretations of the coupled equations (26) and (27) with coefficients 
in front of R.H.S. given by N^- 

Applying the Fierz transformation for the Fock exchange energy term (26) 
we can see that there appear the color-singlet scalar, pseudoscalar, vector and 
axial-vector self-energies. In general we must take into account these self- 
energies mV,V = Us + jbUps + 7fj.U{^ + 7/i75t^ai, with the mean-fields C/j. 
Here we retain only Us,Uy, U^^ in V and suppose that others to be vanished. 
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We shall see this ansatz gives self-consistent solutions for Eq.(26) because of 
axial and reflection symmetries of the Fermi seas under the zero-range approxi- 
mation for the gluon propagator. We furthermore discard the scalar mean-field 
Us and the time component of the vector mean-field for simplicity since 
they are irrelevant for the spin degree of freedom. 

According to the above assumptions and considerations the mean-field V in 
Eq.(24) renders 

V = J375Ua, Ua = UI, (28) 

with VA Ua- Then the diagonal component of the Green function Gii{p) is 
written as 

Gil {p) = [g^^ - 70 A^toGa a] (29) 

with 

Ga\p) = ^ - m + ^ - J513UA, (30) 
Ga^{p) = ^ - m - ^ - "jETsUa, (31) 

where 7573 = 7573 and Ga{p) is the Green function with Ua which is deter- 
mined self-consistently by way of Eq. (26). 

Before constructing the gap function A, we first find the single-particle 
spectrum and their eigenspinors in the absence of A, which is achieved by di- 
agonaUzation of the operator G^^. We have already known four single-particle 
energies e± (positive energies) and —e± (negative energies), which are given 
as 



e± (P) = y p2 + [/2 + ^2 ± 2UA\lm?+pl (32) 

and the eigenspinors <^s, s = ± should satisfy the equation, G'^ieg, v)(t>s = 
0. 

Here we take the following ansatz for A: 

A(p) = ^A,(p)5,(p), 

s=± 

Bs{p) = 7O0-.(P)<^1(P). (33) 

The structure of the gap function (33) is then inspired by a physical con- 
sideration of a quark pair as in the usual BCS theory: we consider here the 
quark pair on each Fermi surface with opposite momenta, p and — p so that 
they result in a linear combination of J'^ = 0^, 1~ (see Fig. 3). ^ 

As is still a matrix in the color-flavor space. Since the antisymmetric nature 
of the fermion self-energy imposes a constraint on the gap function [1], 



CA(p)C-^ = A^(-p). 



(34) 
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Figure 3. Deformed Fermi seas and the 
quark pair on each surface. The top figures 
show those in the absence of A± and the 
middle figures diffusion of the Fermi sur- 
faces in the presence of A±. The bottom 
ones show the quark pairing on the Fermi 
surfaces. 
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Figure 4. Schematic view of the polar- 
angle dependence of the gap functions at 
the Fermi surface, (a) for m = and (b) 
for m ^0. 



An (p) must be a symmetric matrix in the spaces of internal degrees of free- 
dom. Taking into account the property that the most attractive channel of the 
OGE interaction is the color antisymmetric 3 state, it must be in the flavor 
singlet state. 

Thus we can choose the form of the gap function as 

(a,) ^..=e"'^V^A, (35) 

for the two-flavor case (2SC), where a, (5 denote the color indices and i, j the 
flavor indices. Then the quasi-particle spectrum can be obtained by looking for 
poles of the diagonal Green function, Gn: 



p^„^ _1 x/MpIH^F+TOpF for color 1,2 

^'^^^ -\ ^{es(p)-nf for color 3 ^^^^ 
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Note that the quasi-particle energy is independent of color and flavor in this 
case, since we have assumed a singlet pair in flavor and color. 

Gathering all these stuffs to put them in the self-consistent equations, we 
have the coupled gap equations for A^, 

^s'{Ke,) = ^9 j^P '^^^Y.Ts's{K^k.pA):^^^y (37) 
and the equation for Ua, 

within the "contact" interaction, (f' = (/^/A^, (see Eq. (40)), where ^^(p) 
denotes the momentum distribution of the quasi-particles. We find that the 
expression for Ua, Eq. (38), is nothing but the simple sum of the expectation 
value of the spin operator with the weight of the occupation probabiUty of the 
quasi-particles Vg for two colors and the step function for remaining one color 
(cf. (19)). 

Carefully analyzing the structure of the function Ts's in Eq. (37), we can 
easily find that the gap function should have the polar angle (6) dependence 
on the Fermi surface, 

A.(pf , 0) = f-. , R + f] , (39) 

with constants F and R to be determined (see Fig. 4). 

As a characteristic feature, both the gap functions have nodes at poles (9 = 
0, tt) and take the maximal values at the vicinity of equator (9 = tt/2), keeping 
the relation, A_ > A_|_. This feature is very similar to pairing in liquid ^He 
or nuclear matter [17, 18]; actually we can see our pairing function Eq. (39) to 
exhibit an effective P wave nature by a genuine relativistic effect by the Dirac 
spinors. Accordingly the quasi-particle distribution is diffused (see Fig. 3) 

Self-consistent solutions 

Here we demonstrate some numerical results; we replaced the original OGE 
by the "contact" interaction with the cutoff around the Fermi surface in the 
momentum space, 

D^^- ^.gf^'^/A^^ As{p)^As{p)9{S-\es-fi\) (40) 

as in the BCS theory in the weak-coupling limit [19]. 

First we show the magnitude of Ua (Fig. 5). It is seen that the axial-vector 
mean-field (spin polarization) appears above a critical density and becomes 
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larger as baryon number density gets higher. Moreover, the results for different 
values of the quark mass show that spin polarization grows more for the larger 
quark mass. This is because a large quark mass gives rise to much difference 
in the Fermi seas of two different "spin" states, which leads to growth of the 
exchange energy in the axial-vector channel. A slight reduction of Ua arises 
as a result of diffuseness of the Fermi surface due to A^. As seen in Eq. (38), 
Ua can be obtained as addition and cancellation of the contributions by two 
different Fermi seas; the latter term is more momentum dependent than the 
former one and thereby Vg{p) enhances the cancellation term (see Fig. 3). 




Figure 5. Axial- vector mean-field (VA) 
as a function of baryon number density 

psipo = 0.16fm~^). Solid (dashed) lines 
denote VA in the presence (absence) of 
CSC. 



Figure 6. Mean values of the gap func- 
tions, A± and their mass dependence. 



Next we show the gap function as a function of pB (Fig. 6). To see the bulk 
behavior of the gap function, we use the mean-value with respect to the polar 
angle on the Fermi surface. 




The mean values (A±) begin to split with each other at a density where Ua 
becomes finite. We'd like to make a comment here. One may be surprised to 
see their value of O(GeV), coming from our parameter choice. However, what 
we'd Uke to reveal here is not their realistic values but a possibility of color 
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magnetic superconductivity and its qualitative features. More realistic study, 
of course, is needed by carefully checking our approximations, especially the 
contact interaction and the sharp cutoff at the Fermi surface. 

With these figures we can say that FM and CSC barely interfere with each 
other. 

4. Chiral symmetry and magnetism 

We have seen that the quark mass dependence of ferromagnetism should 
be important, while we have treated it as an input parameter. When we con- 
sider the realization of chiral symmetry in QCD, the quark mass should be 
dynamically generated as a result of the vacuum "superconductivity"; qq pairs 
are condensed in the vacuum. We consider here SU{2)l x SU{2)r symme- 
try. Then Lagrangian should be globally invariant under the operation of any 
group element with constant parameters, except the symmetry-breaking term 
stemming from the small current mass, rric. Here we'd like to suggest another 
mechanism leading to FM in quark matter with recourse to chiral symmetry. 
We shall see that FM may be realized without accompanying spin polarization. 

Consider the following parameterization for the combination of the quark 
bilinear fields by introducing the auxiliary fields, p and Oi'. 

-ij^ilj + ij5Ttpij5Ttp = p ex.p{i'y5T ■ 0) . (42) 

Then it resides on the chiral "circle" with "modulus "p and " phase" , any 
point on which is equivalent with each other in the chiral limit, rric = 0, 
and moved to another point by a chiral transformation. We conventionally 
choose a definite point, (vac|/?|vac) = fjr (f-jr' the pion decay constant) and 
(vac|6'j|vac) = 0, for the vacuum, which is flavor singlet and parity eigenstate. 
In the following we shall see that the phase degree of freedom is related to spin 
polarization; that is, the "phase condensation" with a non-vanishing value of 
leads to FM [20]. 

Separating the fields p and Oi into the classical ones and fluctuations around 
them, and discarding any fluctuation the mean-field theory proceeds: 

p-^{p)[>0], e.^iOi). (43) 

Assuming the simplest but nontrivial form of the classical chiral angle such that 
^3(r) = q • r, 9i^2 = 0, we call this set a dual chiral density wave (DCDW) ^ : 

(^■0) = A cos q • r 
{■ipij^Tsip) = Asinq-r. (44) 

It should be obvious that if A vanishes, the phase degree of freedom has to 
become redundant, as seen later. It would be worth mentioning that similar 
configuration has been studied in other contexts [21-23]. Note that the config- 
uration in (44) breaks rotational invariance as well as translational invariance, 
but the latter invarianvce is recovered by an isospin rotation [26]. 
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Taking the Nambu-Jona-Lasinio (NJL) model as a simple but nontrivial 
model [27], we explicitly demonstrate that quark matter becomes unstable for 
a formation of DCDW above a critical density; the NJL model has been origi- 
nally presented to demonstrate a realization of chiral symmetry in the vacuum, 
while recently it been also used as an effective model embodying spontaneous 
breaking of chiral symmetry in terms of quark degree of freedom [28] ^ 

Lnjl = - mc)il^ + G[{'^^f + {i>ii5T^f] (45) 

Using the mean-field approximation (MFA) with the DCDW configuration, we 
introduce a new quark field il)w by the Weinberg transformation [29], 

i)w = exp[i75T3q ■ r/2]^, (46) 

to get the transformed Lagrangian, 

jOmf = i^w[i9-M- 1/275X3^] V'w - GA2, (47) 

with the dynamically generated mass, M = —2GA and = (0, q). This 
procedure embodies translational invariance of the ground state, and shows that 
we essentially consider a "uniform" problem, while we introduced the space- 
dependent mean-fields at the beginning. We briefly summarize in Table 3 the 
relation of the transformed frame to the original one. Note that the transformed 
Lagrangian becomes the same as the familiar form used in discussions of chiral 
symmetry realization within the NJL model, except the isovector and axial- 
vector coupling term 75X3^. We can see the role of the wave vector q is the 
same as AV introduced in the previous sections. 

Table 3. Diagram of the Weinberg transformation. 

ii^ip) ¥=0 ^ i^w^w) = 0) 

(^nsTsV') (i'wi'ysTaiJw) = 

q/2 oc Ve ("AV") 
non-uniform uniform 

The Dirac equation for tpw then renders 

{i^-M- 1/2t375^) = 0. (48) 

We can find a spatially uniform solution for the quark wave function, ij^w = 
uw{p) exp(ip • r), ^ and the energy eigenvalue is given as 

= s]eI + |q| 2/4 ± ^(p • q)2 + M2|qj^, = (M^ + |p|2)V2 (49) 
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for positive energy (valence) quarks with different polarizations. For negative 
energy quarks in the Dirac sea, they have a spectrum symmetric with respect 
to null line because of charge conjugation symmetry in the Lagrangian (47). 
The single-particle spectrum (49) shows again an analogous feature to the ex- 
change splitting between two eigenenergies with different polarizations in the 
presence of q; hereafter, we choose q//z, q = (0, 0, g), q>Q, without loss 
of generality. 

Thus each flavor quark shows the same energy spectrum (49) even in the 
presence of the isospin dependent AV and its form is the same as in Eq. (14). 
However, there is one and important difference from the previous sections; 
we have considered the flavor singlet AV, while we are now considering the 
isovector AV here. The eigenspinor u^ ,-, i = u,d for each flavor satisfies the 
same Dirac equation for a given energy eigenvalue, except the different sign in 
the AV term, so that we have the different form for each flavor; = u^^. 

The mean-value of the spin operator = ( ^ ) ^^en given by 



(73 

-s^ - ^u^^ S n± - (50) 



with (3 = ^/pI + M^ for u quarks. The corresponding value for d quarks is 
also given as s^^ = —sf^. Thus we can see two flavors are oppositely polar- 
ized to each other. Since the integral of sf^ over the Fermi seas should be finite 
for g / for each flavor, the spin polarization of each flavor is finite but has 
opposite direction to each other. Consequently the total spin polarization or the 
flavor singlet AV is always vanished in this case. However, note that this result 
is never conflicted with FM of quark matter by considering the magnetization. 
As we have already noted, the response of the system to the magnetic field 
goes through the magnetization. Taking into account the difference of electric 
charges of two flavors Qi, Qu = +2/3e and Q^, = —l/3e, we can see that 
each flavor coherently contributes to the magnetization, instead. 

Thermodynamic potential 

The thermodynamic potential is given as 



= f^val + ^^vac + MV4G. (51) 

where 6± = 0{i.i — E^)^ M the chemical potential and 7 is the degeneracy 
factor 7 = NfNc- The first term ilvai is the contribution by the valence quarks 
filled up to the chemical potential, while the second term Uvac is the vacuum 
contribution that is formally divergent. We shall see both contributions are in- 
dispensable in our discussion. Once fitotai is properly evaluated, the equations 
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to be solved to determine the optimal values of A and q are 

^^total C^^total 



dA dq 



= 0. (52) 



Since the NJL model is not renormalizable, we need some regularization 
procedure to get a meaningful finite value for the vacuum contribution. Con- 
sider the sum of the negative energy over the Dirac sea, 

J^va. = -7 E / - ^-f, (53) 

where we subtracted an irrelevant constant i^i-cf = — 27/ ■^:^Ep with an ar- 
bitrary reference mass M = M^cf to make the following procedure mathemat- 
ically well-defined. Since the energy spectrum is no more rotation symmetric, 
we cannot apply the usual momentum cut-off regularization (MCOR) scheme 
to regularize $7vac- Instead, we adopt the proper-time regularization (PTR) 
scheme [30]. We think this is a most suitable one for our purpose, since Qvac 
counts the spectrum change under the "external" axial-vector field. It has been 
shown that the vacuum polarization effect under the external electromagnetic 
field can be treated in a gauge invariant way, where the energy spectrum is also 
deformed depending on the field strength [30]. It is also known that the con- 
sequences from the NJL model are almost regularization-scheme independent 
[28], including the PTR scheme. 

Introducing the proper-time variable r, we eventually find 

(54) 

except an irrelevant constant Oref, which is reduced to the standard formula 
[28] in the limit q 0. 

We can easily see, from Eq. (54), that the q degree of freedom becomes su- 
perfluous and theory must become trivial in the limit m ^ 0, which is equiv- 
alent to A in the chiral limit: all the observables must be independent of 
q. This salient feature is consistent with the form of DCDW. The integral with 
respect to the proper time r is not well defined as it is, since it is still divergent 
due to the T contribution. Regularization proceeds by replacing the lower 
bound of the integration range by 1/A^, which corresponds to the momentum 
cut-off in the MCOR scheme. 

For given chemical potential and M and q we can evaluate the valence 
contribution ri^ai using Eq. (49) and write down the general formula analyti- 
cally. Then the thermodynamic potential can be expressed as ilyai = ^vaiil) — 
I^Pvaiiq), where e^aiiq) and Pvai(Q') are the energy density and the quark-number 
density, respectively. They consist of two terms corresponding to the two 



Possibility of color magnetic superconductivity 



17 



Fermi seas with different polarizations: evai('7) = + ^^^^ Pvaiil) = 
Pvaxil) + Ptaii^)- ^® present some examples about the instability of the usual 
NJL ground state with respect to spontaneous generation of DCDW. In the 
present calculation chiral symmetry restoration occurs at the first order in the 
case without DCDW. ^ We can see the NJL ground state becomes unstable at 
the critical chemical potential /ici, and symmetry restoration is delayed until 
fic2 by the presence of DCDW. This dragging effect by DCDW is one of the 
important features. 
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Figure 7. Countour map of the thermo- 
dynamic potential in the dynamical mass 
(Af) - wave number (q) plane. The abso- 
lute minimum is denoted by the cross for 
given density. We have the first order phase 
transitions in this calculation. 



Figure 8. Density dependence of M 
and q, compared with the usual result with 
q = 0. There appear two critical chemical 
potentials; the lower one indicates the in- 
stability of the ground state for formation 
of DCDW, while the higher one restoration 
of chiral symmetry. 



5. Summary and Concluding remarks 

In this talk we have discussed a magnetic aspect of quark matter based on 
QCD. First, we have introduced "ferromagnetism" (FM) in QCD, where the 
Fock exchange interaction plays an important role. Presence of the axial- vector 
mean-field (AV) after the Fierz transformation is essential to give rise to FM, 
in the context of self-consistent framework. As one of the features of the rel- 
ativistic FM, we have seen that the Fermi sea is deformed in the presence of 
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AV; the Fermi sea has a prolate shape for the majority spin particles, while an 
oblate shape for the minority spin particles. 

We have then discussed a possibility of color magnetic superconductivity 
and seen coexistence of FM and CSC is possible. Our ansatz for quark pairing 
shows an effective P- wave pair condensation and gives a polar angle depen- 
dence of the gap function, which looks similar to liquid ^He A - phase. Note 
that this ansatz is never unique for color magnetic superconductivity and other 
types may be also possible [16], where the gap function should show other an- 
gle dependence. In this context recent studies about S = 1 quark pairing may 
be interesting [31]. 

We have briefly discussed a relation of magnetism to chiral symmetry and 
presented an idea, dual chiral density wave (DCDW), which should lead to 
FM. Using, e.g., the NJL model we have demonstrated under what conditions 
the ground state becomes unstable for formation of DCDW. We have found the 
usual ground state becomes surely unstable at the critical density and stays in 
FM between the first and the second critical densities. 

The FM induced by DCDW has many interesting features different from the 
Bloch mechanism. Unfortunately we have not revealed them yet, but it would 
be interesting to examine whether DCDW is possible in the CSC phase. 

The symmetry breaking pattern is summarized as follows: in the condensa- 
tion of the flavor singlet AV, it violates rotation symmetry, 

0(3) ^ 0(2), (55) 

while the DCDW state does flavor symmetry as well as rotation symmetry, 

0{3) X SU{2)v ^ 0{2) xUi^{l). (56) 

The latter situation is similar to the neutral pion condensation in nuclear matter. 

It would be important to figure out the low energy excitation modes (Nambu- 
Goldstone modes) built on the ferromagnetic phase. The spin waves are well 
known in the Heisenberg model [10]. Then , how about our case [32]? 

If quark matter is in the ferromagnetic phase, it may produce the dipolar 
magnetic field by their magnetic moment. Since the total magnetic dipole mo- 
ment Mq should be simply given as Mq = fj,q ■ {Ait/ 3 ■ r^)nq for the quark 
sphere with the quark core radius and the quark number density Uq. Then 
the dipolar magnetic field at the star surface R takes the maximal strength at 
the poles, 

«„.^f(5)%,n,^10.nGl(g)'(^)(^ 



We have not considered the electromagnetic interaction between quarks and 
the induced magnetic field. It would be interesting to see how the situation 
changes when we take it into account; symmetry restoration [33] or mixing 
between magnetic field and gluon field are among them [34]. 



Possibility of color magnetic superconductivity 



19 



Finally we'd like to give a comment about fluctuations. In this talk we have 
completely discarded fluctuations and been only concerned with the mean- 
field. It would be reasonable to study the phase transition, at least qualitatively. 
However, we know some fluctuations or correlations between relevant opera- 
tors should have some effects even before the phase transitions. In particular 
the axial and magnetic susceptibilities in normal quark matter would be inter- 
esting; they might have important consequence,e.g., for quark-quark pairing 
correlation as in ^He superfluidity [17]. 
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Notes 

1. FM does not necessarily accompany SSP in some cases with internal degrees of freedom, as is seen 
in section 4. 

2. There appears no tensor mean field in QCD as a result of chiral symmetry. So we, hereafter, only 
consider AV. 

3. On the contrary, the Fermi sea remains spherical in the non-relativistic case [10]. It would be also 
interesting to compare our results with those given in the different context [12]. 

4. Note that this choice is not unique; actually we are now studying another possibility of quark pair 
between different Fermi surfaces [16]. 

5. Some authors considered similar configuration [24] and called a chiral density wave in analogy 
with spin density wave (SDW) by Overhauser in condensed-matter physics [25]. However, only the scalar 
density oscillates with finite wave number and the pseudo-scalar one is discarded in their ansatz. 

6. We can see that the OGE interaction gives the same form after the Fierz transformation in the zero- 
range Umit. 

7. This feature is very different from refs.[24], where the wave function is no more plane wave. 

8. Note that this is not a unique possibility: we may have the second order phase transitions for other 
parameter sets [28]. 
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